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ABSTRACT 

We investigate the variability of exoplanetary radio emission using stellar magnetic maps and 
3D field extrapolation techniques. We use a sample of hot Jupiter hosting stars, focusing on the 
HD 179949, HD 189733 and r Boo systems. Our results indicate two time-scales over which 
radio emission variability may occur at magnetised hot Jupiters. The first is the synodic period 
of the star-planet system. The origin of variability on this time-scale is the relative motion 
between the planet and the interplanetary plasma that is co-rotating with the host star. The 
second time-scale is the length of the magnetic cycle. Variability on this time-scale is caused 
by evolution of the stellar field. At these systems, the magnitude of planetary radio emission 
is anticorrelated with the angular separation between the subplanetary point and the nearest 
magnetic pole. For the special case of r Boo b, whose orbital period is tidally locked to the 
rotation period of its host star, variability only occurs on the time-scale of the magnetic cycle. 
The lack of radio variability on the synodic period at r Boo b is not predicted by previous 
radio emission models, which do not account for the co-rotation of the interplanetary plasma 
at small distances from the star. 
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1 INTRODUCTION 


Within the solar system, the magnetised planets are all known to be 
sources of intense non-thermal radio emission ( |Zarka||2007t . The 
emitted power is proportional to the power incident on the planet’s 
magnetosphere. Surprisingly, the conversion efficiency between in¬ 
cident and emitted power is roughly constant for every planet. This 
relation is encapsulated by the radiometric Bode’s law (Desch &] 


|Kaiser|[T^ IFarrell et al.|[T9^ |Zarka et al.||200Tl |Zarka||2007> 

with Jupiter having the strongest radio emission of the solar system 
planets. 

In analogy with the solar system, it is thought that magnetised 
exoplanets may emit in the radio wavelengths. However, some stel¬ 
lar systems contain short-period Jupiter-mass planets ( [Bayliss ^ 
|Sackett|2011[[Wright et al.|2012| |. Due to the small star-planet dis¬ 
tances (<0.1 AU), the amount of energy from the solar wind inci¬ 
dent on so called ‘hot Jupiters’ is higher and the magnitude of radio 
emission is expected to be correspondingly larger ( |Stevens|2005| ). 

Numerous authors have modelled exoplanetary radio emission 
with the aim of identifying the most promising targets for future ob- 
servations (e.g. [Lazio et al.|2004[|GrieBmeier et al.|2007b[|Jardine| 
[& Collier Cameron[[2008[ l. Despite repeated attempts, there have 

been no confirmed exoplanetary radio emission detections to date 


E-mail: wcvs@st-andrews.ac.uk 


([Bastian et al.[ 2000| [Ryabov et al.[[2004| [Lazio & Farrell [2007| 

Smith et al.[[2009 

[Lecavelier Des Etangs et al.[[2009[ 

201 1[ Hal-1 

linan et al.|2013| 

Sirothia et al.|2014|). [Bastian et al.[([2000[l offer 


several possible reasons for the lack of detection. The observations 
may have been made at a frequency different to that of the radio 
emission or could simply have lacked the sensitivity required. Al¬ 
ternatively, the radio emission may not have been beaming towards 
Earth. These explanations do not rule out the existence of radio 
emission, just our ability to detect it. The authors also argue that 
these systems could lack a source of keV electrons necessary to 
produce the emission. The lack of a detection would then be due to 
a genuine lack of emission. Lastly, the authors discuss the sporadic 
nature of radio emission. They argue that emissions above their de¬ 
tection threshold may exist but, due to emission variability, it was 
not detectable at the time of their observations. 

With few exceptions ( [Fares et ^[2010[ [Vidotto et al.[[2012[ 
[2015[ ), emission variability has been neglected in radio emission 
models to date. Most assume stellar magnetic fields and winds that 
are both steady and isotropic. However, observations of the solar 
wind (e.g. [Gosling| 1996[ ), and modelling of stellar magnetic fields 
(e.g. [Petit k al.[2008[ l and winds (e.g. [Vidotto et al.[2009[ l show that 
this is not the case. Indeed, the dynamic nature of the solar wind 
is one of the factors, among others, known to affect the intensity 
of radio activity at Jupiter and Saturn ( [Zarka|19^ [Gurnett et al.[ 
[2002[[CTary et al.[2005[ l. 
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In this paper we model the planetary radio emission from 
seven planet hosting systems, whose magnetic fields have been 
mapped by [Fares et ^ ( |2013| ). In particular, we focus on HD 
179949, HD 189733 and r Boo. We incorporate these published 
magnetic maps into our model and consider the time-scales over 
which planetary radio emission variability might be expected. 

The rest of the paper will be structured as follows. In sec¬ 
tion 1^ we outline our radio emission model. The sample of hot 
Jupiter hosting stars is discussed in sectionIn sectionwe dis¬ 
cuss the results of the radio emission modelling for HD 179949b, 
HD 189733b and r Boo b, and the rest of the sample and conclud¬ 
ing remarks are given in section]^ 


2 RADIO EMISSION MODEL 


The model of exoplanetary radio emission employed in this work 
is adapted from the model presented by jJardine & Collier Cameron] 
( |2008| ) and reproduces well the observed ‘radiometric Bode’s law’. 
Radio emissions originate from a population of electrons that have 
been accelerated in the current sheet that forms where the stellar 
and planetary magnetic fields interact. The power of the accelerated 
electrons is given by 

Pe = VcshrunK, (1) 


where Vcs is the volume of the current sheet, hrun is the rate at which 
runaway electrons are generated, and K is the characteristic energy 
to which each electron is acceleratec0 We take the acceleration 
region to be the size of the planetary magnetosphere, r^s, with a 
fixed aspect ratio, a. Its volume is therefore given by • 

The magnetosphere size, determined by pressure 
balance between the wind ram pressure, stellar magnetic pressure 
and planetary magnetic pressure, is given by 




^J^onwindmpvlj-j. -h BliTorb) 


1/6 


rpu 


( 2 ) 


where Bpi is the planetary magnetic field strength, jjiQ is the perme¬ 
ability of free space, n^md is the wind number density, is the 
proton mass, Vgff is the effective wind speed, B^ir^rb) is the stellar 
field strength at the orbit of the planet and is the planetary ra¬ 
dius. The electric field, E, in the acceleration region determines the 
rate of runaway electron generation and has a dependence given by 


P ^effBi,{rorb) 

— oc -, (3) 

^cs 

where is the Dreicer field, Vgff is the effective wind velocity, and 
Hcs is the density within the acceleration region. jJardine & Collier] 
[Cameron] ( ]20081 ) showed that the radiometric Bode’s law for the 
Solar System planets can be reproduced if the density within the 
current sheet is enhanced by a factor, fcomp^ when compared to the 
solar wind density due to compression within this region, such that 
^cs = fcomp^wind- Wc dircct the interested reader to appendix [A| for 
further details of the parameters E, and Ep,. Finally, the radio 
flux density received at Earth, O, is given by 


^ [jardine & Collier Cameron]|2008) state the power in the accelerated elec- 
trons is given by However, the equation the authors give for the 

runaway electron number density, which they denote as Nmn (their Eq. 13), 
is actually the rate at which runaway electrons are generated per unit vol¬ 
ume with units of m“^s“k The correct form is given by our Eq. (^. This 
modification does not affect the scalings, and hence conclusions, drawn by 
these authors. 


O = 


^Id^df 


(4) 


where Prad is the radio emission power, Q is the solid angle into 
which the emission is beamed, d is the distance to the system from 
Earth, and 6f is the bandwidth given by the maximum cyclotron 
frequency, eBpijlnnie ( [Griefimeier et al. [2007a] ). Here, e and are 
the electron charge and mass respectively. Laboratory experiments 
have shown that roughly 10% of the power in the accelerated elec¬ 
trons is converted to radio emission such that Prad - O.IP^. Eor 
this study, we adopt values of K =10KeV, a - 10"^, D=1.6sr and 
fcomp-i^ which are characteristic values for these parameters, and 
Jovian parameters for the exoplanet. 

Two factors control the magnitude and variability of radio 
emissions at each system. The wind density determines how large 
a pool of electrons are available for acceleration within the cur¬ 
rent sheet while the stellar magnetic field determines the amount 
of energy available to accelerate those electrons. Though the wind 
density decreases as a function of orbital distance, roughly as r~^ 
once the wind is terminal, it does not vary significantly around the 
planetary orbit within our model. Instead, it is the stellar magnetic 
field structure that induces variability in the radio emissions. 

A feature of this model is that the radio flux density calculated 
by Eq. 0 is independent of the assumed planetary field strength. 
A stronger field results in a larger planetary magnetosphere (Eq. 
1^). Since the volume of the current sheet is related to the magne¬ 
tosphere size, Vcs ^ ^mp^ ^ larger pool of electrons is available for 
acceleration (Eq. 0)- The increased power is offset by increased 
loses, associated with a higher bandwidth, as the radio signal prop¬ 
agates to Earth (Eq. (|^). Since the frequency of radio emission de¬ 
pends on the planetary magnetic field strength ( ]Zarka]2007] ), which 
are observationally unconstrained, our model does not make a pre¬ 
diction on the frequency of radio flux density at Earth, just its mag¬ 
nitude. 

The power of the accelerated electrons depends on the size 
of the magnetosphere, the effective velocity and the local electric 
field. These are determined by the large-scale flow and magnetic 
field. We discuss these aspects of the model in more detail in the 
following subsections. 


2.1 Stellar magnetic field extrapolation 

In this section, we describe the process by which the 3D magnetic 
field structure of cool stars can be determined using a potential-field 
source surface (PESS) approach ( [Altschuler & Newkirk] 1969| l. The 
magnetic field is assumed to be in a potential state (Vx5 = 0) which 
allows it to be defined in terms of a scalar potential, B = -Vjp- 
Since no magnetic monopoles exist in nature, V • 5 = 0, Laplace’s 
equation arises naturally, V^xj/ - 0. Solving Laplace’s equation in 
terms of spherical harmonics gives us the coefficients of the asso¬ 
ciated Legendre polynomials, P/^, from which the three magnetic 
field components at any point in the corona can be determined: 


N I 

Br = -Yj - (/ + (cos 6) 

/=1 OT=1 

^ ^ d 

l=l m=l 


Be 


N I 


l=l m=l 


(cos 9) 

sine* 


(5) 

( 6 ) 

(7) 
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Figure 1. Examples of the magnetic field geometries of HD 179949 (2009 September), HD 189733 (2008 July), and r Boo (2011 January). For each star, maps 
of the radial magnetic field at the stellar surface (first row) and source surface (second row) are shown. The colour bars indicate the magnetic field strengths 
in Gauss. Polarity inversion lines are shown with dashed lines. White dots on the surface maps indicate the foot points of the field lines that connect with the 
planet around the orbit. The large-scale 3D coronal field structure between the stellar surface and the source surface is also shown (third row). Closed and open 
field lines are coloured in red and blue respectively. Finally, the stellar wind ram pressure, for all latitudes and longitudes, at the orbit of the planet are shown in 
units of 10“^ dyne cm“^ (fourth row). All the stars display a complex field at the photosphere but a dipole dominated field exists at the source surface because 
higher order spherical harmonics decay more quickly with distance above the stellar surface. Three dimensional structure exists in the ram pressure and the 
stellar magnetic pressure (which has a structure similar to that seen in the second row). 


where aim and bim are the amplitudes of the spherical harmonics, / 
is the spherical harmonic degree and m is the order or ‘azimuthal 
number’. 

To determine the values of aim and bim, we apply two boundary 
conditions. Magnetic maps reconstructed using Zeeman-Doppler 
imaging (ZDI) are used to set the first boundary condition at the 
stellar surface, r^. Examples of such surface maps are shown in Fig. 
T1(top row) for HD 179 949 (|Fares et al.|2012> . HD 189733 ( |F^ 
^ al.|2010t and t Boo jFares et al.|2012r ' The second boundary 
condition is set at the source surface, which represents the limit 
of coronal confinement. At the source surface, the magnetic field is 
forced to be purely radial, i.e. Bq - - 0 ([Altschuler & Newkirk 

|1969[|todine et al.|2002| ). Beyond the source surface, the field re¬ 
mains purely radial, decaying as an inverse square law. Physically, 
the source surface represents the location where the gas pressure is 
large enough to blow open closed coronal loops. We adopt a source 
surface distance of = 3.4r^. Since higher order spherical har¬ 
monic modes decay more rapidly with distance away from the stel¬ 


lar surface (Eqs. ID-lz}. the field structure is predominantly dipolar 
at the source surface for these stars. This is true despite the more 
complex fields evident at the stellar surfaces shown in Fig.[2 

To carry out the field extrapolation, we use a modified ver¬ 


sion of the global diffusion code of [van Ballegooijen et al.j l p^^ . 
The PFSS approach has previously been used to study a variety of 


stars (jJardine et al. 

[2002 1 [Gregory et al. [[2006 1 [Lang et aL[[2012| 

[Johnstone et al.|20] 

Examples of the 3D extrapolation for HD 


179949, HD 189733, and r Boo are shown in Fig.[^as well as the 


field geometry at the source surface. 

2.2 Stellar wind 

Equation ^ requires us to estimate two properties of the stellar 
wind. These are the effective wind velocity, i.e. the speed at which 
the interplanetary plasma impinges on the planetary magnetosphere 
from the reference frame of the planet, and the density of the wind. 
Both of these quantities may vary around the orbit of the planet. 
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4 V. See et al. 


For short-period planets, the effective velocity, v^//, has a 
component from the wind velocity, v^ind and a component in the 
azimuthal direction, due to the orbital motion of the planet. 
The effective velocity is calculated by adding these components in 
quadrature. 


Table 1. Properties of the HD 189733, HD 179949 and r Boo systems. 
Listed are the spectral type, stellar mass and radius, planetary orbital 
radius, stellar rotation period, planetary orbital and synodic periods and 
distance from the Sun. All values referenced from |Fares et al.|j2013) 
with the exception of Psyn (calculated) and d (referenced below). 


V 


2 

eff 


2 

wind • 


( 8 ) 


At the small orbital distances where hot Jupiters exist, the in¬ 
terplanetary plasma will be corotating with the star. The azimuthal 
velocity component is therefore given by Vaz - InVorb/Psyn where 
rorb is the orbital radius and Psyn = Pro!-P^b synodic period 

of the planet with respect to the rotation of the star. Prot and Porb 
are the stellar rotation and planetary orbital periods respectively. In 
the rotating frame of the host star, it takes one synodic period for 
a planet to return to its starting position. Consequently, radio emis¬ 
sions are periodic with a time-scale given by the synodic period 
( [Fares et al.|2010| ). 

We calculate using a two step process. The wind speed 
is estimated at the source surface radius first. Solar wind speeds 
are known to correlate with the amount of field line divergence 
( Levine et al.jj 19771 [Wang & Sheeley|ri990j l. [Wang & Sheeleyj 
( 1991| t showed that such a correlation is plausible provided that 
the amount of Alfven wave energy flux is constant within open flux 
tubes. Building on the work of [Wang & Sheeley1 ( |1990j ), p^rge &| 
|Pizzo| ( [2000] ) establish a continuous empirical relation between the 
magnetic expansion factor, fs, and the solar wind velocity at the 
source surface of the Sun, 


410 

Vwindirss) = 267.5 -f ^[ms 

Js 


The magnetic expansion factor is given by 


(9) 


Pss j 


( 10 ) 


where B{rss) is the magnetic field strength at a location on the 
source surface and 5(ro) is the magnetic field strength at the so¬ 
lar surface along the same magnetic field line. Since the so called 
Wang-Sheeley-Arge model is calibrated for the solar wind, we 
modify Eq. replacing r© with r^. B(r^) and B(rss) are deter¬ 
mined in the coronal field extrapolation by eqs. 0-0. 

The second step is to propagate v^mdi^ss) out to Porb- Above the 
source surface, we assume that the wind evolves according to the 
model of |Parker| ( |1958| ). The wind velocity is found by integrating 
the momentum equation for the wind. 


dv^ind ^ e 1 rj.. GM^ 

— = - — {phsT)- pm —(11) 

where p is the mass density, m is the average molecular mass, T 
is the isotropic wind temperature and the other symbols have their 
usual meanings. The wind temperature is chosen to match the ve¬ 
locities at the source surface. The wind velocity at the orbit of the 
planet can then be determined using this temperature. Our approach 
to calculating the wind velocity allows the full 3D structure of the 
stellar magnetic field to be imprinted on the wind velocity. Thus the 
wind ram pressure varies around the planetary orbit in a way that is 
determined by the observed magnetic field geometry. 

Presently, no direct methods of observing the very low density 
winds of low-mass stars exist. |Jardine & Collier Cameron j ( |2008| ) 
use a scaled solar wind density in their model. These authors set 
the wind density to be Nw = 1.7 x at a distance of 

r = 215r^ from the host star. The density at the exoplanet is then 
calculated according to mass conservation and assuming a wind 


Parameter 

HD 189733 

HD 179949 

T Boo 

SpT 

K2V 

F8V 

F7V 

[Mo] 

0.82 

1.21 

1.34 

r* [rol 

0.76 

1.19 

1.42 

rorb [AU] 

0.031 

0.0439 

0.048 

Prot [days] 

12.5 

7.6 

3.3D 

Porb [days] 

2.22 

3.09 

3.31 

Psyn [days] 

2.7 

5.2 

oo 

d [pc] 

19.3^ 

21.V 

15.6“^ 


^Fares et ^ j2013) use the equatorial rotation period of 3 days 
for T Boo. We choose to use a rotation period 3.31 days, which is still 
consistent with the differential rotation exhibited by r Boo, to explore 
the effects of a tidally locked planet. 

^Fares et al.j J2010|, ^Butler et al.j (|2006|, ^GrieBmeier et al.j 
<2007a> 


that is terminal and radial, i.e. the wind density falls as r~^. How- 
ever, since we use a radially varying wind in this work, the density 
falls as in order to conserve mass. Additionally, we scale 

the wind density by a factor, f^ag, where fmag is the average foot 
point strength of the field lines that intersect the exoplanet’s orbit 
normalised to the average Solar surface field (IG). Physically, fmag 
accounts for the denser winds of more active stars ( jMestel & Spruitj 
|1987| ). The density at the orbit of the planet is therefore given by 


P^windiP) 


10^ f 

Jmag 


Vwindjr = 215r^) 

^windif) 



( 12 ) 


In the last row of Fig.[^ we show the ram pressure of the wind on 
a shell at the orbital radius of the planet. The influence of the mag¬ 
netic field structure can clearly be seen in the ram pressure struc¬ 
ture. 


3 STELLAR SYSTEM PROPERTIES 

In this paper, we use the sample of seven hot Jupiter hosting stars 
presented by [Fares et al.| ( |201^ for which a ZDI map exists. In par¬ 
ticular we focus on HD 179949, HD 189733 and r Boo. Properties 
for these three systems are shown in Table[^with details for the full 
sample available in [Fares et al.[ ( [2013'] ). 

We choose these three systems for further discussion because 
they are often cited as promising targets for detectable exoplanetary 
radio emission ([Lazio et al.[2004[[Griefimeier et al.[2007b[[Jardine[ 

[& Collier Cameron[2008[ ) and have been extensively studied in the 
literature. In addition, these stars allow a comparison of two criti¬ 
cal factors that affect exoplanetary radio emission. The first is the 
spectral type of the host star which affects the nature of the stel¬ 
lar wind. HD 179949, an F8 dwarf, and HD 189733, a K2 dwarf, 
allow us to make this comparison. The second is the relative angu¬ 
lar velocity of the planet’s orbital motion and the stellar rotational 
motion which determines whether the planet moves relative to the 
stellar magnetic field. We can evaluate the effect of this factor with 
HD 179949b and r Boo b. Both planets have host stars with similar 
spectral types but the latter system is known to be tidally locked 
whilst the former is not. 

HD 189733b and r Boo b are known to orbit in, or close to. 
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Figure 2. Properties of the HD 179949 system during the 2007 June (left) and 2009 September (right) epochs. Top: Magnetic field maps of the source surface 
displayed in the same format as Fig.^ Bottom: Predicted planetary radio flux (red) and stellar magnetic field strength along the planetary orbit (blue) as a 
function of the stellar longitude of the subplanetary point. Large drops in radio emission are evident as the planet crosses the polarity inversion line. 


the stellar equatorial plane. Making use of the Rossiter-McLaughlin 
effect ( |Rossiter|1924||McLaughlin|1924| ), [TTraud et al.| ( |2009| ) de¬ 
termine the spin-orbit m isalignment angle f or HD 189733b to be 
0.85°Q 28 ^- For 'r Boo b, Brogi et al. |2012|) find that the normal 
to the orbital plane is not significantly misaligned with the stel¬ 
lar rotation axis. For the purposes of this study, we will assume 
that both planets lie exactly in the stellar equatorial plane. For HD 
179949b, [Brogi et ar](|2014|) deter mine that the orbital inclination 
is 67.7° + 4.3°. [Fares et al.| ( |2012| ) determine that the stellar incli¬ 
nation is roughly 60° using ZDI. Inclinations determined from ZDI 
can have errors of 10° or greater (e.g. [Petit et ar]|2008| ). We will 
assume a spin-orbit misalignment of 0° for HD 179949 since this 
is consistent with the stellar and orbital inclinations when consid¬ 
ering the errors. [Winn et ^ ( [2010j ) note that hot Jupiter systems 
with hotter host stars {Tgff > 6250K) are more likely to have a 
large spin-orbit misalignment. Within our sample, HD 179949 and 
T Boo are the hottest stars, with effective temperatures of 6168K 
and 6387K respectively ( [Fares et al.|2013| l, and are the most likely 
to have large spin-orbit misalignments according to the trend noted 
by [Winn et al.| ( |20l'0| ). Nevertheless, we will proceed by assuming 
the planetary orbits are equatorial for the reasons given above. For 
the remaining five systems, the spin-orbit misalignments are un¬ 
known since the planets are not transiting. For the purposes of this 
study, we assume that they also orbit within the stellar equatorial 
plane. 


In terms of the geometry of the system, it is also important to 
consider the fact that the radio emission is anisotropic. The emis¬ 
sions are beamed into a cone which is aligned with the local plan¬ 
etary magnetic field ( [Zarka|2007] l. If the magnetic dipole axis and 
planetary rotation axis are aligned, the radio emissions should be 
beamed in the direction of Earth at all times, if the geometry is 
favourable, or not at all. However, if there is a misalignment be¬ 
tween the planetary magnetic dipole and rotation axes, then emis¬ 
sions, as observed from Earth, may be modulated on the planetary 
rotation period as the emission cone swings in and out of view. 
Given that we have no information about the planetary field align¬ 


ment, we will assume the most favourable geometry, in which the 
radio emissions are always beamed towards Earth. 

Eor each star in our sample, magnetic maps are recon¬ 
structed from spectropolarimetric observations. The observations 
were made using one of two high-resolution spectropolarimeters, 
either ESPaDOnS at the 3.6-m Canada-Erance-Hawaii Telescope 
(CEHT) on Mauna Kea or NARVAL at the 2-m Telescope Bernard 
Lyot (TBL) in Pic du Midi (Prance). Two magnetic maps are avail¬ 
able each for HD 179949 ( [Fares et al.|20T2l ) and HD 189733 ( [F^ 
jet al.|2010]). For r B oo, we use seven ZDI maps spanning four years 
( [Donati et al. [2008 [ [Fares et al.|2009||2013j ). ZDI maps for the re¬ 
maining four stars are available in Fares et al. j pOT^ . In Fig. [2 
we show examples of the magnetic field structure for HD 179949, 
HD 189733, and r Boo. In each case, the radial component of the 
magnetic field at the stellar surface & source surface is shown, as 
well as the 3D coronal field structure between these two bound¬ 
aries. Though complex structure is evident at the surface of each of 
the stars, it is predominantly the dipole component that survives at 
the source surface. 


4 RESULTS 
4.1 HD 179949 

Figure[^shows magnetic maps at the source surface of HD 179949 
(top row). These maps are obtained using the extrapolation process 
described in section [TT] using magnetic maps of the stellar surface 
obtained during 2007 June and 2009 September ( [Fares et al.[2012[ |. 
The colours represent the magnetic field strength in Gauss and the 
dashed lines represent polarity inversion lines, i.e. locations where 
Br = 0. To calculate the expected planetary radio emission, we 
require the local stellar magnetic field strength in the vicinity of 
HD 179949b as it orbits over one synodic period. This is obtained 
by considering the stellar magnetic field strengths in the equatorial 
plane of the source surface maps and calculating the magnetic field 
strength at the distance of the planet according to an inverse square 
law decay. The local field strength has been plotted in the bottom 
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Figure 3. Properties of the HD 189733 system during the 2007 June (left) and 2008 July (right) epochs. Data displayed in the same 


format as Fig.|^ 


row of Fig. in blue as a function of the stellar longitude of the 
planet. The corresponding predicted radio flux received at Earth, 
over one synodic period, is also plotted in red. Figure has been 
formatted such that the stellar longitudes of the top and bottom 
rows are aligned. 

Over each synodic period (5.2 days), HD 179949b exhibits 
highly variable radio emission that correlates strongly with the lo¬ 
cal stellar held strength. The variability can be mostly attributed 
to the orbital motion through the stellar magnetic held. Comparing 
the source surface maps with the predicted radio flux curves, neg¬ 
ligible emission occurs as the planet crosses the polarity inversion 
lines. Conversely, the highest emission occurs when the planet is 
over the magnetic poles, i.e. when the held strength is strongest. 
This is true during both epochs, between which, the stellar held has 
evolved. 

As well as planetary radio emission, several authors have sug¬ 
gested that magnetic interactions may also induce enhanced chro¬ 
mospheric activity near the stellar surface ( [Shkolnik et al.|[2008l 
|Lanza|200^ . [Fares et al.| ( [?012| ) searched for this type of interac¬ 
tion at HD 179949 flnding hints of activity enhancement modulated 
at the synodic period of the system. Given that enhanced chromo¬ 
spheric activity and planetary radio emission are both thought to 
occur as a result of magnetic interactions between a planet and host 
star, it is unsurprising that both would be modulated on the synodic 
period. 


4.2 HD 189733 

In Fig.[^ we plot the same information as Fig. but for the HD 
189733 system. The ZDI maps used for the held extrapolations 
were presented by [Fares et ^ ( [2010[ ). Comparing Figs. and 
we see that HD 179949b and HD 189773b both display strong vari¬ 
ability, modulated on the synodic period of their respective systems 
during a given epoch. Given that both planets orbit stars where the 
dipole component of the magnetic held is tilted, the similarity is 
to be expected. The main difference between these systems is the 
magnitude of the radio flux density, smaller by roughly an order of 
magnitude at HD 179949b. This can be attributed to the lower pho- 


tospheric magnetic held strength, the larger star-planet distance, 
and the greater distance from Earth for HD 179949. 

Using the model presented by [Griefimeier et al.] ( [2007b[ ), [Fares[ 
[et al.[ ( [20T0] ) have also predicted the radio emission expected from 
HD 189733b. The shape of their radio flux curve as a function of 
subplanetary point is in agreement with ours though the magnitude 
of their radio emission is higher, reaching the hundreds of mJy, due 
to a higher conversion efflciency from the Poynting flux to radio 
emission. 

A number of authors have searched for radio emission from 
HD 189733b. While none produced positive detections, 3cr upper 
limits of 2.1 mJy, 2mJy, 81mJy, and 160//Jy were found at fre¬ 
quencies of 150MHz ([Lecavelier Des Etan gs et al.[2011| ), 244MHz 
([Leca velier Des Etan gs et al.[[2009|l, 307-347MHz ([Smith et al.[ 

[2009 j ), and 614MHz ( [Lecavelier Des Etangs et~^[2009[ t respec- 

tively. Recalling that our model does not predict emission of a spe¬ 
cific frequency, we can compare these limits to our predicted radio 
flux densities. Assuming our predictions are correct, it is unsurpris¬ 
ing that [Smith et al.[ ( [200^ were not able to detect any radio emis¬ 
sions given their upper limit of 81mJy. However, at the other three 
frequency bands, the upper limits are either comparable, or sig¬ 
nificantly lower, than our predictions. [Lecavelier Des Etangs et ^ 
( [2009[ ) discuss several reasons for the lack of detection and con¬ 
clude that the most likely is because the emission was at a lower 
frequency. 


4.3 T Boo 

Similarly to Figs. and we plot magnetic maps of the source 
surface and the predicted planetary radio emission for the r Boo 
system in Fig.[^ We use seven ZDI maps, with epochs spread over 
three years, for the held extrapolations of r Boo ([Donati et al.[2008| 
[Fares et al.[200^[20T^ . 

Unlike the previous systems, the orbital motion of r Boo b is 
tidally locked to its host star ( [Donati et al.[2008[ l. Consequently, the 
subplanetary point is fixed at one stellar longitude (180° in Fig.|^. 
Since there is no relative motion between planet and host star, ra¬ 
dio emission variability from planetary motions through the stellar 
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Figure 4. Properties of the r Boo system over seven epochs. Data is displayed in the same format as Figurej^ The predicted radio flux displayed is that which 
could be expected if r Boo b’s subplanetary point sampled every stellar latitude. However it remains fixed at a single stellar longitude, 180°, due to tidal 
locking. Black vertical lines are included to indicate this longitude. The expected radio flux, at any given epoch, is given by the intersection of the red and 
black lines. Figure[^shows how this flux varies with time. 


field no longer exists. The only source of variability, in this model, 
is from dynamo driven evolution of the stellar field between epochs. 
It is worth highlighting that this behaviour would not be seen from 
models where the stellar rotational velocity is neglected, such that 
the Vaz component of Vg/f is set to the Keplerian velocity of the 
planet. If Vaz is non-zero, there must be some relative motion be¬ 
tween the planet and the surround interplanetary plasma and mag¬ 
netic field. 

Given the lack of relative motion between host star and planet. 


we would not expect to see the entirety of the predicted radio curves 
in Fig.|^ Rather, these curves represent the expected radio flux if r 
Boo b were not locked to its host star. In reality, we would only ex¬ 
pect to see the radio flux indicated at 180° during each epoch, which 
has been highlighted by black vertical lines. In Fig.[^ we show the 
predicted radio flux from r Boo b as a function of time by collecting 
into one plot the radio fluxes, at 180° for each epoch, from Fig.|^ A 
different symbol is used for each epoch to allow easier comparison 
between Figs. s-lzl This plot demonstrates that variability occurs 
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Figure 5. Predicted radio flux of r Boo b as a function of time. The con¬ 
necting lines are visual aids only and should not be interpreted as a form 
of interpolation. Each epoch is plotted with a different symbol to allow eas¬ 
ier comparison with Figs. and [7] In contrast to HD 189733, where the 
principle source of variability can be attributed to the orbital motion of the 
planet, variations in radio flux are caused by dynamo drive In evolution of 
the stellar magnetic held. 



Figure 6. Location of the negative pole of the dipole component of r Boo’s 
magnetic held during each epoch. Lines connecting contiguous epochs are 
visual aids only and should not be interpreted as a form of interpolation. For 
each epoch, the symbol used is the same as that used in Fig.[^ The negative 
poles during 2009 June and 2010 January epochs are very close (latitude 
~8° and longitude ~180°) and only appear to be a single data point in this 
plot. 



Angular separation (°) 


Figure 7. Radio flux as a function of angular separation between the sub¬ 
planetary point and the closest magnetic pole of r Boo during each epoch. 
Symbol shapes correspond to those used in Figs. [^and[^ Filled symbols are 
used if the nearest pole is negative (the points plotted in Fig.|^ and unfilled 
symbols are used if the nearest pole is positve (antipodal to those in Fig.[^. 
A trend of decreasing radio flux can clearly be seen for increasing angular 
separation from the nearest stellar magnetic pole. 


Table 2. Maximum predicted planetary radio fluxes for the rest of the 
systems presented by [Fares et al.| ( |2013) . 


HD No. 

102195 

130322 

46375 

73256 

^max [mJy] 

0.062 

6.0 X lO-'* 

0.018 

0.19 

Psyn [days] 

6.2 

18.2 

3.3 

3.1 


over months/years, i.e. the time-scale of the field evolution. The 
predicted radio fluxes are of the same order of magnitude as those 
predicted by jVidotto et al.| ( [2012| ) despite the differing models and 
assumptions used. 

The radio flux variations in Fig.j^occur due to evolution of the 
stellar held. The magnitude of the radio flux is highest when r Boo 
b is over a magnetic pole, e.g. 2008 January, and weakest when it is 
over a magnetic cusp, e.g. 2008 June. Figure shows the location 
of the negative pole of the dipole component {l-l mode) during 
each epoch (red markers) and the subplanetary point (blue marker). 
Figure shows the predicted radio flux as a function of angular 
separation between r Boo b’s subplanetary point and the nearest 
magnetic pole during each epoch (filled symbols if the nearest pole 
is negative and unfilled if it is positive). A clear trend of decreasing 
radio flux for increasing angular separation to the nearest magnetic 
pole is apparent. This result highlights the importance of knowing 
the three dimensional structure of the stellar magnetic held and not 
just the structure in the orbital plane of the exoplanet. It also ar¬ 
gues for further ZDI observations of r Boo which may allow us to 
determine any periodic behaviour in its large-scale magnetic held. 
If the evolution of the dipole component could be predicted, some 
estimate of the radio flux in the future could be made. 


4.4 Rest of sample 

From the sample of [Fares et al.| ( |2013] ), four further stars had de¬ 
tectable magnetic fields with one ZDI map existing for each. None 
of these planets are tidally locked to their host star. If any of these 
planets are radio emitters, all should show variability as a result of 
their orbital motions. Tableshows the predicted maximum plan¬ 
etary radio flux received at Earth per synodic period and the syn¬ 
odic period for each of these systems. None of these planets are 
particularly strong radio emitters. [Vidotto et al.[ ( [2015[ ) have also 
calculated the radio fluxes emitted from these exoplanets. With the 
exception of HD 130322, they And fluxes which are roughly com¬ 
parable (within a factor of 1-4). The reasons we And such low radio 
fluxes are the same as those given to explain HD 179949b’s lower 
radio flux density compared to HD 189733b’s, namely low pho- 
tospheric magnetic field strengths, large star-planet distances, and 
large distances to Earth. The combination of these factors means 
that these systems are not ideal targets for future planetary radio 
observations. 


5 DISCUSSION & CONCLUSION 

We have presented a planetary radio emission model that incorpo¬ 
rates a realistic geometry for the large-scale stellar magnetic held. 
Additionally, our model estimates the effective wind velocity, v^//, 
more accurately than existing models in the literature. An accurate 
estimate should include a wind component, v^ind, that accelerates 
radially, and an azimuthal component, that accounts for the ex¬ 
oplanetary orbital motion relative to the stellar rotation. [Stevens[ 
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( |2005| ) assumes Ve/f = v^md = 400kms"^ but do not include an az¬ 
imuthal component. This is likely to be an overestimate of and 
an underestimate of VQ^. |Jardine & Collier Camer^ ( |2008| ) assume 
Veff - Vaz for plaucts orbiting close to their host star. These authors 
account for the stellar rotation but do not include a wind compo¬ 
nent. [GneBmeieret^ ([2007^ incorporates a radially accelerating 
wind component and an azimuthal component. However, they do 
not account for the stellar rotation and so overestimate values for 
Vaz and hence Veff. These examples demonstrate the need for care¬ 
ful consideration of the stellar system dynamics when calculating 
the effective velocity parameter. 

Focusing on the HD 179949, HD 189733 and r Boo planetary 
systems, we looked at the time-sc ales on which exoplanetary ra¬ 
dio emissions may vary. We find that planetary radio emissions are 
strongly dependent on local stellar magnetic field strength along the 
exoplanetary orbit. Therefore, exoplanetary radio emissions could 
be used as a probe of the stellar magnetic field geometry. 

In general, the height of magnetic loops on the surface of 
cool stars are comparable to their foot point separation. Therefore, 
higher order multipoles will decay more rapidly with height above 
the stellar surface. These higher order modes will contribute to the 
variability in the radio emissions of exoplanets that are orbiting 
close to the stellar surface of their host star. However, exoplanets at 
larger orbital radii, such as those in our sample, will predominantly 
experience the dipole component of the stellar field. Emission max¬ 
ima exist when the planet is above a pole while a minima would be 
expected inside a magnetic cusp. In general, the magnitude of radio 
emission is anticorrelated with the angular separation between the 
subplanetary point and the closest magnetic pole. 

We find two time-scales for variability. The shorter time-scale 
is the synodic period of the planetary system and is caused by the 
orbital motion of the planet. From the reference frame of the planet, 
the local stellar magnetic field strength is constantly varying caus¬ 
ing corresponding variations in the radio emission, typically on the 
order of days. This type of behaviour, which is exhibited by the HD 
179949 and HD 189733 systems, has implications for future ob¬ 
servations. It is similar to the behaviour described by [Llama et al.| 
( |2013| ) with respect to observing early ingresses for transiting exo¬ 
planets. In both cases, unfavorable wind conditions may lead to a 
lack of detection and highlights the importance of considering the 
field/wind conditions when planning future observations. 

The longer time-scale is most apparent in systems where the 
orbital and stellar rotation periods are similar in length. This re¬ 
sults in a very long synodic period. In the most extreme case, when 
the stellar rotation period and orbital period are equal, there is no 
relative motion between planet and the stellar atmosphere in the az¬ 
imuthal direction. For such systems, e.g. r Boo, variability only oc¬ 
curs as a result of stellar magnetic field evolution. Previous models 
which estimate the azimuthal component of the effective velocity 
to be the Keplerian speed of the planet do not predict this type of 
behaviour for tidally locked systems. Indeed, by definition, the ef¬ 
fective velocity in such systems cannot have a non-zero azimuthal 
component. One would expect variability at locked systems to oc¬ 
cur over a characteristic time given by the magnetic cycle period. 
Non tidally locked systems will display variability on this time- 
scale as well but the time-scale for variability as a result of orbital 
motions will be much shorter. 

Future exoplanetary radio emission detection attempts should 
focus on exoplanets orbiting close to stars with with strong mag¬ 
netic fields and which are close to Earth. Additionally, the exo¬ 
planet would ideally have a large orbital velocity in the rotating 
frame of its host star. This would result in a higher azimuthal com- 
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ponent of the effective wind velocity and hence higher Poynting 
flux. 
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APPENDIX A: RUNAWAY ELECTRON GENERATION 
AND THE ACCELERATION REGION 

The component of the electric field that is parallel to the magnetic 
field in the acceleration region is given by 

£ =(Al) 
F(l/4) 

where H is the magnetic Reynolds number and the remaining sym¬ 
bols are defined in the main text. This is the full version of Eq. 
[Jardine & Collier Cameron| ( [2008) give a full derivation and discus¬ 
sion of this expression. In a plasma, the motions of the electrons, 
under the action of an electric field, are opposed by a drag force 
caused by coulomb collisions. Por electrons with high enough ve¬ 
locities, the coulomb drag becomes negligible and the electrons can 
freely run out of the distribution. The exact behaviour is determined 
by the strength of the electric field. Below a critical field strength, 
known as the Dreicer field ( |Dreicer|1959| t, only the electrons in the 
high velocity tail of the distribution runs away. However, at super- 
Dreicer strengths, the entire distribution is able to do so. The Dre¬ 
icer field is given by 

= 18 X (A2) 

where the rics and T are the electron density and temperature re¬ 
spectively inside the current sheet. Due to compression of the 
plasma within the current sheet, the density does not necessarily 
take on their wind values. Both the compressed density and tem¬ 
perature in the current sheet are unknown. We find that an increase 
in the density by a factor of 15 and a temperature of IMK repro¬ 
duces the emitted radio power at the solar system planets well and 
adopt these values. The rate at which runaway electrons is given by 
[Kruskal & Bernstem| ( |1964| ): 

= 0.35n,,vJ(E/ED) = 2.6 X lO-^nlj-^/^fiE/Eo) (A3) 

where is the electron collision frequency and f^ElEo) is given 
by 
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